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Abstract 

Data on pp — y nnw for beam momenta 600-1940 MeV/c are presented. The strongest 
channel is fo(1500)n from the initial pp state 1 Sq. Together with 777r°7r° data, the 3n data 
determine the branching ratio RR[/ (1500) -)• nn] /BR [/ (1500) -> vr 7r ] = 0.42 ± 0.09. 
They are consistent with a dominant contribution from an / = 0, C = +1 J = 
resonance observed earlier in the 777r°7r° data; from the combined 777r°7r° and nnn data, its 
mass is M = 2320 ±15 MeV and its width T = 230 ± 35 MeV. 

As part of a study of pp annihilation in flight to neutral final states, we have earlier presented 
results on pp — y 777r°7r° [1], with statistics of typically 70,000 events per momentum. We have 
at the same time collected statistics of up to 192 events per momentum on pp — y r\r\r\ in the 67 
channel. Despite the limited statistics, some useful conclusions may be drawn. 

In 7]7T 7c° data, the /2(1270)r/ final state is dominant; the contribution from / (l°00)7r is 
small at all beam momenta: up to 3.2% with errors of ~ 0.6% at each momentum. However, 
the branching ratio of /o(1500) to 7777 is a factor ~ 45 larger than for / 2 (1270); consequently 
/o(1500) makes the larger contribution to 3rj data. These data therefore provide a valuable check 
on some features of the most recent amplitude analysis of pp — y ?77r 7r ; this has recently been 
analysed in a combined fit with data on pp — y 7r°7r°, 7777, rjrj' and 7r~7r + [2]. A 0~ resonance found 
at ~ 2285 MeV in that analysis is predicted to be the main feature of the /o(1500)r7 channel 
here. The present 3r] data allow an improvement in the determination of its mass and width. 

The data were taken at LEAR using the Crystal Barrel detector. Experimental details have 
been given in Ref. [1], and here we consider only features which involve identification of the 377 
final state. In order to isolate this rare channel in 67 data, tight selection criteria are needed 
to eliminate backgrounds from other channels. In processing data, we first demand exactly 
six photon showers, each confined to a block of 3 x 3 adjacent Csl crystals of the detector. 
Extra energy deposits produced by Compton scattering out of the primary showers into nearby 
crystals (so-called 'split-offs') appear in ~ 50% of events. All attempts to recover such events 
lead to backgrounds higher by typically a factor 3-6. We therefore reject events containing such 
additional energy deposits. Events where two photons from one ir° merge into a single shower 
are also discarded. 

A least squares kinematic fit is first required to pp —y 67, satisfying energy-momentum 
conservation with confidence level CL > 10%. Then 7C fits are made to all 3-particle final 
states involving 7r°, 77 or 7/. We reject events which fit 37r°, 7/7r°7r°, rjrjn , 7]'ii tt or 777/71-° with 
CL > 10~ 4 . As further rejection against the prolific 37r°, 777r°7r° and 777771"° channels, events are 
rejected if they fit 7r°7r°77 or 7r°r/77 with CL > 10 -4 . 
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Channel Background (%) 


rjrjn 


0.1 


777r 7r 7r 


0.1 


rjf]7l 7Y 


1.1 


curjrj, oj — > 7r°7 


0.9 


curjn , uj — > 7r°7 


0.1 



Table 1: Background levels in 67 data from competing channels at 1800 MeV/c. 



A Monte Carlo study using 20,000 generated events for all competing channels shows that 
the main sources of potential background are 777771"°, f]ii 7T 7i , i]i]7r 7r , and urjir or urji] where 
oj 7r°7. In the last four cases, two photons or one go undetected. The background levels at 
1800 MeV/c from competing channels are estimated from this Monte Carlo simulation and are 
summarised in Table 1. The principal backgrounds arise from ur/i] and r]r]iT 7T . The situation 
is similar at other momenta and the total surviving background is 2.3% within errors at all 
momenta. Kinematic fits to the 777777 hypothesis are consistent with this estimate and show 
only a few scattered events with M(77) outside the 77 — > 77 peak. 



Momentum 


377 


e 


(7(377) 


(MeV/c) 


events 






600 


9 


0.130 


3.0 ±0.7 


900 


68 


0.126 


4.3 ±0.4 


1050 


72 


0.126 


5.0 ±0.5 


1200 


157 


0.124 


6.9 ±0.5 


1350 


103 


0.121 


6.4 ±0.5 


1525 


99 


0.117 


8.3 ±0.6 


1642 


142 


0.115 


9.8 ±0.6 


1800 


192 


0.113 


9.5 ±0.5 


1940 


163 


0.112 


5.7 ±0.4 



Table 2: Numbers of 377 events after background subtraction and reconstruction efficiency e; the 
last column shows the weighted mean of the integrated cross section for 3i] (corrected for all 77 
decays) from 67 and IO7 data. Errors include statistics and the uncertainty in normalisation of 
cross sections at individual nomenta from Ref. [3] 

Columns 2 and 3 of Table 2 show the number of surviving 377 events and also the reconstruction 
efficiency e. This efficiency is estimated by generating > 50, 000 Monte Carlo events at every 
momentum. These events are subjected to identical selection procedures to data and are used 
for the maximum likelihood fit described below. Cross sections are derived from the number 
of observed events, corrected by the efficiencies of Table 2 and using the number of incident 
antiprotons and the target length of 4.4 cm. A correction is applied for an observed dependence 
of reconstruction efficiency as a function of beam rate; this is described in full in Ref. [3]. 

As a check, the IO7 data have also been examined for the final state 7777371-°, where one 
77 — > 3n°. Events are selected by demanding exactly 10 photon showers and a kinematic fit to 
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Figure 1: The rj — > 37r° peak in pp — » r]r]3Tc data at 900 MeV/c. 



r]r]37r° with CL > 10%. Events fitting 5n° or 77471° with CL > 0.1% are rejected. Unfortunately, 
the background level under the rj — > 3ir° peak is substantial (20% of the signal at the highest 
momenta, rising to 80% at the lowest). It arises mostly from 7]1]3tc° and cannot be reduced. It is 
illustrated at 900 MeV/c in Fig. 1. This background level is too high to allow a physics analysis. 
However, these events may be used, after subtracting background, to check the evaluation of 
integrated cross sections. Fig. 2 compares cross sections for pp — > 3i] from 67 events (circles 
joined by the full curve) and IO7 events (squares joined by the dashed curve). There is agreement 
between these two determinations within the errors. The last column of Table 2 shows the 
weighted mean of the integrated cross sections from 67 and IO7 data. A detail is that statistics 
at 600 MeV/c are much lower than at other momenta; this accounts for the small number of 
events, despite a similar cross section to that at 900 MeV/c. 

Dalitz plots from 67 data and projections on to M(r)T)) are shown at 8 beam momenta from 
1940 to 900 MeV/c in Figs. 3 and 4. Processes which ned to be considered are: 



pp -> /o(1500)77 (1) 

<"7 (2) 

/ 2 (1270)77 (3) 

/o(1370)77 (4) 

/2(1525)77 (5) 



Of these, the last two turn out to be negligible. The /o(1500) is conspicuous at the higher 
momenta. The ^(1270) is obscured in mass projections of Figs. 3 and 4 by reflections from 
the stronger /o(1500). The fit shows that it makes a small but significant contribution. The 
arj channel makes a strong contribution to r]ii ii data, so its presence in 3r] is necessary; it 
accounts well for the small uniform component over the Dalitz plots of Figs. 3 and 4. Here a is 
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Figure 2: Comparison of integrated cross sections from 67 data (circles and full curve) with that 
from IO7 data (triangles and dashed curve). 



a shorthand for the / (400 - 1200) of the Particle Data Group (PDG) [4]. It is fitted with the 
parametrisation of Zou and Bugg [5]. 

The data have been fitted by the maximum likelihood method. Present data are closely 
related to those for pp — y r]7i TT , where statistics are much higher. Our strategy is to take 
magnitudes and phases of all partial waves for /2(1270)?7, fo{1500)r) and ar\ from our analysis of 
77"7r 7r° [2], allowing one free parameter to determine the branching ratio ri27o = BR[f 2(1270) — > 
r)T)]/BR[f 2(1270) — > 7r°7r°], a second to determine the branching ratio ri 50 o = BR[f (1500) — > 
rjrj] /BR[ / (1500) — y 7r°7r°] and a third to determine the branching ratio r a = BR[a — y rjr]\/BR[a — y 
7r°7r ]. The contribution from / 2 (1270) is small, as shown below in Fig. 6(b) by the chain curve. 
The value of r 12 7o can therefore be set to the PDG value. The absolute value of the integrated 
cross section for 3rj is sensitive to r 150 o and r a , where therefore need to be fitted. 



Momentum 


-/o(1500)*7 


-/ 2 (1270)77 


+/o(1370)t? 


+^(1525)t? 


900 


-27.6 


-4.4 


1.1 


1.2 


1050 


-25.9 


-10.8 


1.8 


0.9 


1200 


-44.9 


-1.0 


1.0 


0.1 


1350 


-56.6 


-7.1 


2.2 


0.8 


1525 


-15.0 


-19.8 


0.8 


1.7 


1642 


-21.8 


-41.6 


4.4 


2.0 


1800 


-95.2 


-20.6 


0.3 


0.5 


1942 


-50.6 


-11.5 


4.8 


1.1 



Table 3: Changes in log likelihood (a) removing / (1500)r/ from the fit, (b) removing / 2 (1270)?7, 
(c) adding / (1370)t7 and (d) adding ^(1525)^. 
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Figure 3: (a)-(d) Dalitz plots of data at beam momenta of 1940, 1800, 1642 and 1525 MeV/c, 
(e)-(h) fitted Dalitz plots, (i)-(l) projections on to M(rjrj), (m)-(p) production angular distri- 
bution for events in a band 120 MeV wide centred on / (1500). Histograms show the fit. 
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Columns 2 and 3 of Table 3 show changes in log likelihood when /o(1500) or /2(1270) is 
removed from the fit at individual momenta. Our definition is such that a change of log likelihood 
of 0.5 corresponds to a one standard deviation change for one degree of freedom, so the observed 
changes indicate the presence of both channels. Dropping / 2 (1270) from the overall fit to all 
momenta, log likelihood is worse by 90.4. 

Histograms on Figs. 3 and 4 display results of the fit to all momenta simultaneously. The first 
column of Dalitz plots shows data and the second column the fits. Every event is plotted in three 
7777 combinations. The low statistics cause substantial fluctuations for data and are responsible 
for some apparent disagreement with the fit; the latter is much smoother than data since Monte 
Carlo statistics of over 5000 events per momentum are used in the maximum likelihood fit. 
In fact, the x 2 between data and fit on the Dalitz plot averages 1.1 per point, so there is no 
discrepancy. 

The partial wave amplitude for production of the fa between particles 1 and 2 will be used 
as an illustration of the way amplitudes are parametrised. It takes the form 

f= M ,_ s l_ iMr B^)BMZ(p,k). (6) 

Here G is a complex coupling constant, the denominator refers to the /2(1270), s is mass squared 
for the 771772 combination, and Z is a relativistic Zemach tensor given explicitly in Ref. [1]. The 
cross section is obtained from the coherent sum of amplitudes for the three 7777 combinations. 
The value of G for 3i] data is related to that for i]7i°7i° by G(3i]) = ^ ryxj^k^j k r] G{rj'K Q 'K Q ) / v^3. 



The factor l/y/3 allows for the three identical 7777 pairs; the ratio ^Jk n /k v , involving momenta 
k of 7r and 77 in the decay / 2 — > nn and 7777, allows for phase space for those decays. Standard 
Blatt-Weisskopf centrifugal barrier factors B L with radius 0.8 fm are used to parametrise the 
dependence on the centre of mass momentum p with which the resonance is produced, with 
orbital angular momentum L, and also on its decay momentum k in the resonance rest frame. 

The complex coupling constant G assigns a phase to each partial wave. This phase originates 
from the initial state pp interaction and from rescattering in the final state. The approximation 
we adopt is that the phase is the same for decays / 2 (1270) — > r\r\ and / 2 (1270) — > tt°tt° and 
likewise for /o(1500) — > rjr] and 7r°7r°. This is in the spirit of the isobar model, where it is 
assumed that the final resonant state is reached after any number of intermediate rescatterings 
and then decays into the 7777 and 7t°n channels with the same phase. Relative magnitudes of 
/ 2 (1270) contributions in different partial waves are taken from the fit to r]-K ir data of Ref. 
[2]; likewise for / (1500) and a contributions. However, rescattering is possible in r]-K°iT to 
additional final states a 2 (1320)7r° and ao(980)7r°; therefore we allow an overall phase difference 
of / 2 (1270)t7 relative to /o(1500)?7 and at] at all momenta between i]n 7r° data and 377. The fit to 
377 is therefore related to that to r/^ir by 4 parameters: two relative phases between / 2 (1270)?7, 
/o(1500)ry and err] and two branching ratios r 150 o and r a , common to all momenta. Because the 
contribution from / 2 (1270) is small, results concerning / (1500)r7 have little sensitivity to the 
detailed partial wave decomposition of / 2 (1270)ry. This is because correlations between the two 
channels arise only where bands cross on the Dalitz plots and involve rather few events. 

The prediction from 777r°7r° data [2] is tha the dominant contribution to 377 will be from L = 
decays of a 0~ resonance at 2285 MeV. A small L = 2 component is predicted, roughly 15% of 
the intensity of L = 0. Contributions with L — 1 and 3 are found in Ref. [2] to be very small. 
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Figure 5: The predicted Dalitz plot at 1800 MeV/c for 3 P 1 ->■ / O (1500)?7. 



For / (1500)?7 final states, the angular distributions against centre of mass production angle r 
are shown in Figs. 3 and 4. Fortuitously, the L = and L = 2 amplitudes are nearly orthogonal, 
with the result that interferences between them are small. Hence angular distributions in Fig. 
3 and 4 are close to isotropic. The discrepancy between data and fit in Figs. 4(o) is associated 
with the very small number of events in the /o(1500) mass region. 

A point with interesting consequences is that the three / (1500) bands cross at the centre of 
the Dalitz plot near a beam momentum of 1800 MeV/c. The strong constructive interference 
observed in Fig. 3(b) between the three bands requires that the process goes largely through 
a single pp partial wave with even L. Amplitudes for production with L = 1 cancel at the 
intersection point of the three bands and their sum changes sign about this point. To see 
this, consider / (1500) production from the initial state 3 Pi. For initial helicity m — 0, the 
Clebsch-Gordan coefficient for coupling to pp is zero. For m — 1, the amplitude is proportional 
to psinrexp(i0), where r is the polar angle for production of the resonance and is the 
associated azimuthal angle. This amplitude is therefore proportional to px + ipy, where px,Y 
are transverse momentum components of the resonance. For the coherent superposition of the 
three /o(1500)?y combinations, the resultant at the intersection of the three bands is zero, by 
momentum conservation. The same result extends to L = 3 for / (1500)r] and to production 
of /2(1525) with L — 1 and 3, though the amplitude then contains additional factors for the 
f 2 decay. Fig. 5 shows the predicted Dalitz plot for 3 P 1 ->■ / (1500)r] at 1800 MeV/c. This 
distinctive pattern is absent from the data of Fig. 3. The amplitude analysis confirms that any 
L = 1 or 3 processes producing / (1500) or /2(1525) at 1642, 1800 and 1940 MeV/c are absent 
or very weak (summed cross sections < 10% of L = 0). 

The essential physics conclusions of the amplitude analysis are illustrated in Fig. 6. Fig. 6(a) 
shows a fit (with details given below) to the integrated 3r] cross section averaged between 67 and 
IO7 data. The peak requires, as predicted, a dominant contribution from a J p = 0~ resonance 
in the / (1500)r/ channel; its intensity is shown by the full curve in Fig. 6(b). It appears at a 
slightly higher mass than predicted: 2328 ± 16 MeV compared with M = 2285 ± 20 MeV of 
Ref. [2]. The shift is compatible with the combined errors of the two analyses. The width is 
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Figure 6: (a) The fit to the mean integrated cross section from 67 and IO7 data; (b) the con- 
tributions from 7/(2320) ->■ [/ (1500)r/] L=0 (full curve), arj (dashed), r/ 2 (2267) ->■ [f (1500)rj] L=2 
(dotted), and f 2 (1270)r] (chain curve). 



discussed below. 

In addition, a err] contribution is required, as shown by the dashed curve. This interferes 
destructively with [fo(l500)r]] L=0 . Thirdly, there is a small contribution, shown dotted, from 
r/ 2 (2267) — > [/o(1500)?y]i = 2. The f 2 (1270)r] channel makes a small contribution < 1.3 fib for all 
masses (chain curve). 

The f 2 rj contribution contains too many partial waves and is too small to allow a useful 
determination of the branching ratio r^ro between 7777 and 7r°7r°. This parameter is therefore 
fixed at the PDG value, 1.35 x 10~ 2 [4]. The corresponding ratio ri 50 o for / (1500) may be fitted 
over a range of values, because of interference of this channel with arj. The width of the 0~ 
resonance at 2330 MeV correlates strongly with ri 50 o- Small r values require a narrow width 
for the resonance, so as to reproduce the peak cross section. Amsler et al. [6] determine from 
pp — > r/rj-K and 3n° data at rest a ratio ri 50 o = 0.47 ± 0.21; using this value, the width of the 
0" resonance optimises at T = 240 MeV, distinctly smaller than the value 325 ± 30 MeV found 
from t/ttV data [2]. The WA102 collaboration [9] finds ri 500 = 0.54 ± 0.09, compatible with 
Ref. [6] and with our determination below. Abele et al. [6] find r 150 o = 0.23 ± 0.04, which 
requires T = 154 MeV for the 0~ resonance. This narrow width gives a significantly poorer log 
likelihood in fitting 3rj data and is hard to reconcile with r]7i°7i° data. 
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We find that the fit to 777r°7r° data will adjust towards these new values of M and T with only 
a minor change in log likelihood and in other partial waves. The optimum fit to combined 3i] 
and r]n°TT° data gives M = 2320 ± 15 MeV, T = 230 ± 35 MeV and 

r 1500 = 0.39 ± 0.09. (7) 

These values are used in the curves of Fig. 6. Intensities derived from each component are 
shown in Fig. 6(b) including interferences between all three 7777 channels and after integrating 
over the Dalitz plot. In comparing with corresponding intensities for the 777r7r channel, shown in 
Fig. 2 of Ref. [2], one must take note of two points. Firstly, strong constructive interferences 
between the three 7777 channels in 377 enhance the peaks of Fig. 6(b), for example by 60% for 
[/o (1500) 77] l = o. Secondly, Fig. 6(b) uses a width of 230 MeV for the resonance, while Fig. 2 of 
Ref. [3] uses 325 MeV; this increases the peak height in Fig. 6(b) by a further factor 1.9. 

The <t?7 intensity in Fig. 6(b) is ~ 20% of that for r?7r 7r°. If the a is dominantly non-strange, 
one expects from the composition of the 77, namely 

\rj >=0.8\(uu + dd)/V2> -0.Q\ss >, 

a ratio between 7777 and 7r°7r° decay of (0.8) 4 < k v /k w >; here k V;7T are the momenta in the a 
rest frame in 7777 and nn decays respectively. Averaging over phase space for the a, one finds 
< k v /kit >~ 0.6, and hence a predicted intensity for arj in 3rj of 25% of that in r]-ir n , close to 
the fit. 

From present data, we are unable to demonstrate the presence of the J p = 2~ resonance 
decaying to [fo(^00)r]] L=2 . Deviations from isotropy in angular distributions of Figs. 3 and 
4 are too small to demonstate its presence, though angular distributions are compatible with 
prediction within the large errors. We remark that the 2~ resonance is observed clearly in data 
on pp — > i]'7r°7r°, where it makes a dominant contribution in decays to /2(1270)t/ [9]. 

The fourth and fifth columns of Table 3 show changes in log likelihood when either f (1370)i] 
or ^(1525) is added to the fit with L = 0. Statistically one expects an improvement of 1 
for two extra fitted parameters. The evidence for either / (1370) or ^(1525) being present 
is insignificant. Let a be the decay angle with respect to the beam in the rest fram of the 
resonance. The decay angular distribution (3 cos 2 a — 1) for f 2(1525) produced from the initial 
state 1 _D 2 with L = is distrintive and the amplitude would interfere with f (1500)r]. The L = 1 
production amplitude likewise gives a distinctive contribution to the Dalitz plot, similar to that 
of Fig. 5. The mean fitted value of its cross section is 15% of f (1500)rj but some or all of this 
is undoubtedly statistical noise in the fit. 

From earlier analysis of pp — > 777771"° and 37r° at rest, the ratio BR[fo(1370) — > rji]] / BR[f o(1370) — > 
= 0.056 ± 0.04 [6]. Despite the sizeable error, there was no doubt of the presence of some 
/o(1370) — > r/r] signal in the 777777° data [7], considerably higher than that for / 2 (1270). From 
the small / (1370)t7 signal in 377, we therefore deduce that the signal observed in t]tt -k must be 
almost entirely f2{1270)i] with very little f (1370)i]. This is a useful check on the analysis of 
777r°7r°, where a small / (1370) signal is hard to identify in the presence of a dominant / 2 (1270)7] 
contribution. 

In summary, despite low statistics, four results emerge from this analysis. 

• The 377 and r]n°7T° data together require an I = 0, C = +1, J PC = 0~ + resonance with 
M = 2320 ± 10 MeV, T = 230 ± 35 MeV, decaying to / (1500)?7. 
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• The ratio of intensities in 3i] and r]Tr n gives ri 50 o = 0.39 ± 0.09. 

• The ratio of arj intensities in 3rj and r]iT iT is close to that expected for non-strange 
composition of a = / (400 - 1200). 

• The absence of / (1370) — > r\r\ in 3r\ data is a useful check that it makes very little 
contribution to r]ii 7i data. 
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